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Abstract 

We present a statistical analysis of time-resolved spontaneous emission decay curves from en- 
sembles of emitters, such as semiconductor quantum dots, with the aim to interpret ubiquitous 
non-single-exponential decay. Contrary to what is widely assumed, the density of excited emitters 
and the intensity in an emission decay curve are not proportional, but the density is a time-integral 
of the intensity. The integral relation is crucial to correctly interpret non-single-exponential decay. 
We derive the proper normalization for both a discrete, and a continuous distribution of rates, 
where every decay component is multiplied with its radiative decay rate. A central result of our 
paper is the derivation of the emission decay curve in case that both radiative and non-radiative 
decays are independently distributed. In this case, the well-known emission quantum efficiency can 
not be expressed by a single number anymore, but it is also distributed. We derive a practical de- 
scription of non-single-exponential emission decay curves in terms of a single distribution of decay 
rates; the resulting distribution is identified as the distribution of total decay rates weighted with 
the radiative rates. We apply our analysis to recent examples of colloidal quantum dot emission 
in suspensions and in photonic crystals, and we find that this important class of emitters is well 
described by a log-normal distribution of decay rates with a narrow and a broad distribution, re- 
spectively. Finally, we briefly discuss the Kohlrausch stretched-exponential model, and find that 
its normalization is ill-defined for emitters with a realistic quantum efficiency of less than 100 %. 
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FIG. 1: Schematic of the relation between decay of an excited state X* to the ground state X and 
experimental observable parameters. The density of emitters in the excited state is equal to c(t) 
which can be probed by transient absorption. The emitted light intensity as a function of time 
f(t) is recorded in luminescence decay measurements. In photothermal measurements the released 
heat (g — f){t) after photoexcitation is detected. g{t) describes the total decay, i.e., the sum of the 
radiative and the non-radiative decay. 

I. INTRODUCTION 

Understanding the decay dynamics of excited states in emitters such as semiconductor 
quantum dots is of key importance for getting insight in many physical, chemical and biolog- 
ical processes. For example, in biophysics the influence of Forster resonance energy transfer 
on the decay dynamics of donor molecules is studied to quantify molecular dynamics^. 
In cavity quantum electrodynamics, modification of the density of optical modes (DOS) is 
quantified by measuring the decay dynamics of light sources. According to Fermi's 'Golden 
Rule' the radiative decay rate is proportional to the DOS at the location of the emitter 3 . 
Nanocrystalline quantum dots^^, atoms 6 ' 7 and dye molecules^ are used as light sources 
in a wide variety of systems. Examples of such systems are many different kinds of pho- 
tonic materials, including metallic and dielectric mirrors^ 6 -^ 1 ^, cavities^, metallic filmsAi^, 
twcA 3 -^-, and three-dimensional photonic crystals. 
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Figure [T] shows how observable parameters are related to the decay of an excited state 
X* to the ground state X. In photoluminescence lifetime measurements the decay of 
the number of excited emitters is probed by recording a photoluminescence decay curve 
(f(t)). The number of excited emitters c(t) can be probed directly by transient absorp- 
tion measurements'" and non-radiative decay [g — f)(t) can be recorded with pho- 
tothermal techniques" (see Fig. [I]). g(t) is here defined as the total intensity, i.e., the 
sum of the radiative and non-radiative processes. In this paper we discuss photolumines- 
cence lifetime measurements, which are generally recorded by time-correlated-single-photon- 
countingi. The decay curve f(t) consists of a histogram of the distribution of arrival times 
of single photons after many excitation-detection cycles 1 . The histogram is modelled with 
a decay-function from which the decay time of the process is deduced. 

In the simplest case when the system is characterized by a single decay rate T, the decay 
curve is described by a single-exponential function. However, in many cases the decay is 
much more complex and strongly differs from single-exponential decay^" 1 ^'^ 1 ^'^ 1 ^. This 
usually means that the decay is characterized by a distribution of rates instead of a single 
rats'. For example, ensembles of quantum dots in photonic crystals experience the spatial 
and orientational variations of the projected LDOS explaining the non-single-exponential 
character of the decay'. It is a general problem to describe such relaxation processes which 
do not follow a simple single-exponential decay. Sometimes double- and triple-exponential 
models are justified on the basis of prior knowledge of the emitters^. However, in many cases 
no particular multi-exponential model can be anticipated on the basis of physical knowledge 
of the system studied and a decision is made on basis of quality-of-fit. 

Besides multi-exponential models, the stretched-exponential model or Kohlrausch 
function^ is frequently applied. The stretched-exponential function has been applied to 
model diffusion processes', dielectric relaxation', capacitor discharge^, optical Kerr effect 
experiments- 1 and luminescence decay" 1 ^. The physical origin of the apparent stretched- 
exponential decay in many processes remains a source of intense debate" 1 ^. 

Surprisingly, in spite of the rich variety of examples where non-single-exponential decay 
appears, there is no profound analysis of the models available in the literature. Therefore, we 
present in this paper a statistical analysis of time-resolved spontaneous emission decay curves 
from ensembles of emitters with the aim to interpret ubiquitous non-single-exponential decay. 
Contrary to what is widely assumed, the density of excited emitters c(t) and the intensity in 
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an emission decay curve (f(t) or g(t)) are not proportional, but the density is a time-integral 
of the intensity. The integral relation is crucial to correctly interpret non-single-exponential 
decay. We derive the proper normalization for both a discrete, and a continuous distribution 
of rates, where every decay component is multiplied with its radiative decay rate. A central 
result of our paper is the derivation of the emission decay curve f(t) in case that both 
radiative and non-radiative decays are independently distributed. In this most general case, 
the well-known emission quantum efficiency is also distributed. Distributed radiative decay 
is encountered in photonic media 26 , while distributed non-radiative decay has been reported 
for colloidal quantum dots 3 -^ 1 ^ and powders doped with rare-earths^. We derive a practical 
description of non-single-exponential emission decay curves in terms of the distribution of 
total decay rates weighted with the radiative rates. Analyzing decay curves in terms of 
distributions of decay rates has the advantage that information on physically interpretable 
rates is readily available, as opposed to previously reported analysis in terms of lifetimes. We 
apply our analysis to recent examples of colloidal quantum dot emission in suspensions and 
in photonic crystals. We find excellent agreement with a log-normal distribution of decay 
rates for such quantum dots. In the final Section, we discuss the Kohlrausch stretched- 
exponential model, and find that its normalization is ill-defined for emitters with a realistic 
quantum efficiency of less than 100 %. 



II. DECAY MODELS 



A. Relation between the concentration of emitters and the decay curve 

A decay curve is the probability density of emission which is therefore modelled with a 
so-called probability density function^. This function tends to zero in the limit t — > oo. 
The decay of the fraction of excited emitters ^Q- at time t' is described with a reliability 
function or cumulative distribution function ^1 — f^^ 39 - Here c(0) is the concentration of 
excited emitters at t' = 0. The reliability function tends to one in the limit i — ► oo and 
to zero in the limit t — > 0. The fraction of excited emitters and the decay curve, i.e., the 
reliability function and the probability density function^ 9 -, are related as follows: 
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FIG. 2: (color online) Plot of a non-single-exponential decay of the fraction (black solid curve, 
left axis) and the corresponding total intensity decay curve g(t) (red dashed curve, right axis). The 
curves that describe the fraction of excited emitters and the corresponding intensity decay curve 
are strongly different. In this example, ^jl is the Kohlrausch stretched-exponential decay of the 
fraction (Eq. \19\ black solid curve) and g{t) the corresponding decay curve (Eq. [20l red dashed 
curve). We have taken /?=0.5 and T s t r = 1. 

Physically this equation means that the decrease of the concentration of excited emitters 
at time t' is equal to the integral of all previous decay events. Or equivalently: the total 
intensity g(t) is proportional to the time-derivative of the fraction of excited emitters. As 
an illustration, Fig. [2] shows a non-single-exponential decay function simultaneously with 
the corresponding decay curve. It is clear that both curves are strongly different. In many 
reports, however, the distinction between the reliability function and the probability density 
function is neglected: the intensity of the decay curve g(t) is taken to be directly proportional 
to the fraction of excited emitters This proportionality only holds for single-exponential 
decay and not for non-single-exponential decay, which has important consequences for the 
interpretation of non-single-exponential decay. 
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FIG. 3: (color online) Luminescence decay curve of emission from a dilute suspension of CdSe 
quantum dots (open dots, left axis). Data were collected at the red side of the emission maximum 
of the suspension, at A = 650 ± 5 nm. Single-exponential modelling (red dashed curve, right 
axis) yields a decay time of 39.0±2.8ns and a \r °f 1-12. The average photon arrival time < t >, 
calculated with Eq. is 39.1 ns. 

B. Single- exponential decay 

In this section, we will illustrate some concepts with the well-known single-exponential 
model. We will also indicate which features of single-exponential decay will break down in 
the general case of non-single-exponential decay. It is well known that in case of first-order 
kinetics the rate of decrease of the concentration is constant in time: 

dc(t') , ,. 

where T is the decay rate of the process. As a consequence, the concentration c(f') decreases 
single-exponentially in time: 

g = « p(-«') (3) 

The mathematical expression for the luminescence decay curve can be obtained by inserting 
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Eq. [3] into Eq. [H where T is identified with the total decay rate T to t, resulting in: 



g(t) = T rad exp(-T tot t) + T nrad exp(-T tot t) 



(4) 



where T rad is the radiative decay rate, T nrad is the nonradiative decay rate and T tot is the 
total decay rate with Y tot = T rad + T nrad . In a luminescence decay measurement the recorded 
signal is proportional to the first term of g(t) only which is f{t): 



and therefore a single-exponential luminescence decay process is modelled with Eq. The 
pre-exponential factor a is usually taken as adjustable parameter, and it is related to several 
experimental parameters, i.e., the number of excitation-emission cycles in the experiment, 
the photon-collection efficiency and the concentration of the emitter. Henceforth a will be 
omitted in our analysis. A comparison between Eqs. [5] and [3] shows that in the case of pure 
single-exponential decay neglect of the distinction between the reliability function (Eq. [3]) 
and the probability density function (Eq. [5]) has no important consequences, since both the 
fraction and the decay curve are single-exponential. As Fig. [2] shows, this neglect breaks 
down in the case of non-single-exponential decay. 

Figure [3] shows a luminescence decay curve of a dilute suspension of CdSe quantum dots 
in chloroform at a wavelength of A = 650 ± 5 nmM, with the number of counts on the 
ordinate and the time on the abscissa. Clearly, the data agree well with single-exponential 
decay as indicated by the quality-of-fit Xr °f 1.12, close to the ideal value of 1. This means 
that all individual quantum dots that emit light in this particular wavelength-range do so 
with the same rate of ns -1 . It appears that the rate of emission strongly depends on the 
emission frequency and that it is determined by the properties of the bulk semiconductor 
crystal. 

Since f(t) as given by Eq. [5] is a probability density function, the probability of emission 
in a certain time-interval can be deduced by integration. The total probability for emission 
at all times between t = and t — > oo is given by 



which is equal to the luminescence quantum efficiency. The luminescence quantum efficiency 
is defined as the probability of emission after excitation^. The correct recovery of this result 
in Eq. [6] shows that Eq. [5] is properly normalized. 



f(t) = aT rad exp(-T tot t) 



(5) 
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The average arrival time of the emitted photons or the average decay time can be calcu- 
lated by taking the first moment of Eq. [5j 

<(>=T = £-/(«)** = J_ (7) 

Only in the case of single-exponential decay the average decay time < t > is equal to the 
inverse of the total decay rate T to t- The average arrival time for the data in Fig. [3] was 
< t >= 39.1 ns, very close to the value of 39.0 ±2. 8ns obtained from single-exponential 
modelling, which further confirms the single-exponential character of the decay of quantum 
dots in suspension. 

C. Discrete distribution of decay rates 

In contrast to the example shown in Fig. [3], there are many cases in which decay curves 
cannot be modelled with a single-exponential function. As an example, Fig. H] shows a 
strongly non-single-exponential decay curve of spontaneous emission from CdSe quantum 
dots in an inverse opal photonic crystal 15 ' 26 . If a non-single-exponential decay curve is 
modelled with a sum of single-exponentials, the decay curve has the following form: 

1 11 

f(t) = -7-7 C i T rad,i eXp(-T totti t) (8) 

H U J i=i 

where n is the number of different emitters (or alternatively the number of different envi- 
ronments of single emitters^), Cj is the concentration of emitters that has a radiative decay 
rate Y ra d,i-, and c(0) is the concentration of excited emitters at t = 0, i.e., the sum of all con- 
centrations Q. When the different fractions (or environments) are distributed in a particular 
way, a distribution function p(T tot ) may be used. Such a function describes the distribution 
or concentration of the emitters over the emission decay rates at time t = 0. The fraction 
of emitters with a total decay rate T to t,i is equal to 

Ci 1 {c(T tot:i _i) + c(T to t,i+i)) 



c(0) c(0) 



(9) 

, p(r t ot)dr tot 

A JT t ot,i-l 

= p(r tot}i )Ar tot 

where p{T tot)i ) expresses the distribution of the various components i over the rates T tot i 
and has units of inverse rate s. AT tot is the separation between the various components i in 
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FIG. 4: (color online) Luminescence decay curve of emission from CdSe quantum dots in a titania 
inverse opal photonic crystal (dots, left axis). The lattice parameter of the titania inverse opal 
was 340 nm and the emission wavelength A = 595nm. (a) A log-normal distribution of rates (Eq. 
[T71 and [TSJ re d dashed curve, right axis) models the data extremely well {Xr — 1-17). Xhe T m y is 
91.7 /xs _1 (j^-y = 10.9 ns) and the width of the distribution Ar is 0.57 ns -1 . (b) In contrast, 
a Kohlrausch stretched-exponential model (red dashed curve, right axis) does not fit the data 
(Xr = 60.7). The stretched-exponential curve corresponds to r<, ir =96.2 /xs -1 (pl-=10.4 ns), an 
average decay time < t > of 31.1 ns, and a /3-value of 0.42. 

the sum. The decay curve now has the following mathematical form: 

n 

f(t) = AT tot p(T totji )r rad:i exp(-r toM t) (10) 
1=1 

It is important to note that in Eq. [TfJ every component in the sum is correctly normalized 
since every component is multiplied with its radiative decay rate T ra ^i. 
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D. Continuous distribution of decay rates 



For infinitesimal values of AT to t, Eq. [10] can be written as an integral: 

roo 

f(t) = T rad (T tot )p(T tot )exp(-T tot t)dT tot (11) 
Jo 

In the case of single-exponential decay the distribution function is strongly peaked around a 
central T tot - value, i.e., the distribution function is a Dirac delta function. Inserting a Dirac 
delta function into Eq. [TT] recovers Eq. 

roo 

fit) = / T rad , 8 (Ttot - Ttot') exp(-T tot t)dT tot 
Jo 

= T rad , exp(-T tot d) (12) 

This result confirms that the generalization to Eq. [11] is correct since it yields the correctly 
normalized single-exponential functions. 

In Eq. [TTJit is tacitly assumed that for every T tot there is one T rac i: the function T rac i(T tot ) 
relates each T tot to exactly one T rad . In general both T tot and T rad vary independently, and 
Eq. [11] is generalized to 

/(*) (13) 



rco rTtot 

/ / dT rad pr tot (T rad ) T rad 
Jo Jo 



p(T to t) exp(-T tot t)dT 



tot 



where pr tot (r ra d) is the normalized distribution of T rad at constant T tot . For every T tot 
the integration is performed over all radiative rates; a distribution of T rad is taken into 
account for every Ttot- Eq. [13] is the most general expression of a luminescence decay curve 
and a central result of our paper. From this equation every decay curve with a particular 
distribution of rates can be recovered. An example described by Eq. [13] is an ensemble of 
quantum dots in a photonic crystal. In photonic crystals the local density of optical states 
(LDOS) varies with the location in the crystal and the distribution of dipole orientations 
of the emitters^. Therefore, an ensemble of emitters with a certain frequency emit light 
with a distribution of radiative rates T rad . In addition, when an ensemble of emitters has 
a distributed T tot and a single radiative rate T rad , i.e., pv tot (T rad ) is a delta-function, then 
Eq. [T3] reduces to Eq. [TT] Even though the non-radiative rates may still be distributed, 
Eq. [TT] suffices to describe the decay curve since for every T tot there is only one T rad . Such a 
situation appears, for example, with powders doped with rare earth ions^ and with polymer 
films doped with quantum dots^ 1 ^. 
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FIG. 5: Log-normal distribution of T. This distribution was modelled to the data of Fig. 0] (curve 
a, quantum dots in photonic crystal) and Fig. [3] (curve b, quantum dots in a diluted suspension), 
with T m f and Ar as adjustable parameters. For (a) r m j is 91.7 us^ 1 (j^— = 10.9 ns) and the 
width of the distribution Ar was 0.57 ns -1 and for (b) F m f was 25.8 us^ 1 (j^j = 38.8 ns) and 
the width of the distribution Ar was 0.079 ns^ 1 . 

Interestingly, an ensemble of emitters with a distribution of rates T tot is not completely 
characterized by a single value of the quantum efficiency (as opposed to Eq. [6] for single- 
exponential decay). In such an ensemble, the quantum efficiency is distributed, since each 
T tot is associated with a distribution of radiative rates T ra( i. The average quantum efficiency 
< QE > can be calculated by integrating Eq. [13] for all times: 



Most often, detailed information on the relation between T tot and T ra( i is not available. 
Then, modelling directly with a distribution of decay rates is applied^i'^^'^L^. This ap- 
proach has a major advantage over modelling with a stretched-exponential function, where 




(14) 
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it is complicated to deduce the distribution of decay rates (see below). A function of the 
following form is used to model the non-single-exponential decay curve: 

/(*)=/ a(T tot )ex P (-r tot t)dT tot (15) 
Jo 

In Eq. [15] the various components are not separately normalized as in Eq. [131 Modelling 
with Eq. [15] boils down to using an infinite series of single-exponentials which are expressed 
with only a few free parameters. The form of the distribution can usually not be predicted 
and a decision is made on basis of quality-of-fit. While a good fit does not prove that the 
chosen distribution is unique, it does extract direct physical information from the non-single- 
exponential decay on an ensemble of emitters and their environment 

It is widely assumed that cr(T) is equal to the distribution of total ra.tesr±-^££&k£i. A com- 
parison with Eq. [T3l shows that this is not true and reveals that cr(T) contains information 
about both the radiative and non-radiative rates: 

r^tot 

a(T tot ) = p(T tot ) pYtoX^rad) T rad dT rad (16) 
Jo 

Thus cr(r) is the distribution of total decay rates weighted by the radiative rates. This 
conclusion demonstrates the practical use of Eq. [13] the equation allows us to completely 
interpret the distribution of rates found by modelling with Eq. [15] Such a complete inter- 
pretation has not been reported before. 



E. Log-normal distribution of decay rates 

Distribution functions that can be used for <y(T) are (sums of) normal, Lorentzian, and 
log-normal distribution functions. In Fig. S^a) the luminescence decay curve of quantum 
dots is successfully modelled with Eq. [T5] with a log-normal distribution of the rate T 



cr(r) = Aexp 



lnT -lnT 
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(17) 



where A is the normalization constant, T m f is the most frequent rate constant (see Fig. ED- 
7 is related to the width of the distribution: 

Ar = 2T mf sinh( 7 ) (18) 

where Ar is equal to the width of the distribution at -. The most frequent rate constant T m f 
and 7 are adjustable parameters, only one extra adjustable parameter compared to a single- 
exponential model. Clearly, this model (Eq. [T5land [TTj) describes our non-single-exponential 



13 



experimental data extremely well. The xl was 1-17, T m / was 91.7 us' 1 = 10.9 ns) and 
the width of the distribution Ar was 0.57 ns~ l . In addition to T m f and Ar, an average 
decay rate can be deduced from the log- normal distribution in Fig. [5] . However, this average 
is biased since the various components are weighted with their quantum efficiency, as shown 
in Eq. M 

Modelling with a log-normal distribution of decay rates yields direct and clear physical 
parameters, for instance the shape and width of the decay rate distribution. The log-normal 
function is plotted in Fig. [5] (curve a). The broad distribution of rates demonstrates the 
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strongly non-single-exponential character of the decay curve. In Ref. 
relate the width of this broad distribution to the spatial and orientational variations of the 
LDOS in inverse-opal photonic crystals. 

The log-normal model was also modelled to the decay curve from quantum dots in sus- 
pension (Fig. [3]). The distribution is plotted in Fig. [5] (curve b). T m f was 25.8 fis^ 1 
{jt-j — 38.8 ns), close to the lifetime deduced from the single-exponential modelling of 39.0 
±2. 8ns. The narrow width of the distribution Ar of 0.079 ns -1 is in agreement with the 
single-exponential character of the decay curve. 

F. Stretched-exponential decay 

Besides the multi-exponential models discussed in Sections C-E, the Kohlrausch 
stretched-exponential decay mode&^ 9 is widely applied to model non-single-exponential de- 
cay curves. The fraction of excited emitters, i.e., the reliability function, of the Kohlrausch 
stretched-exponential model is equal to: 

= exp(-(r sir ty) (19) 

where (3 is the stretch parameter, which varies between and 1, and T stT the total decay 
rate in case of stretched-exponential decay. The stretch parameter (3 qualitatively expresses 
the underlying distribution of rates: a small (3 means that the distribution of rates is broad 
and P close to 1 implies a narrow distribution. The recovery of the distribution of rates in 
case of stretched-exponential decay is mathematically complicated^ and only feasible for 
specific p>f&&&as.. 

The decay curve corresponding to a Kohlrausch stretched-exponential decay of the frac- 



14 



tion can be deduced using Eq. [19] and Eq. [TJ and results in: 

g(t) = ^(T str tfexp(-(T str tf) (20) 

The normalization of Eq. [20] can, in analogy with Eq. [6] be deduced by integration for all 
times between t = and t — > oo, which yields 1. Therefore, an important consequence is 
that Eq. [20] is correctly normalized only for emitters with a quantum yield of 1 (r ra ^ = T tot 
and f(t) = g(t)). It is not clear how normalization should be done in realistic cases with 
quantum yield < 100%. To the best of our knowledge, this problem has been overlooked in 
the literature. 

The main advantage of modelling with a Kohlrausch stretched-exponential function is 
that the average decay time < t > can readily be calculated. The average decay time is 
equal to22: 

where T is the mathematical Gamma-function. For the single-exponential limit of (3 —> 1 
Eq. [2T] reduces to Eq. [7] Note again that in this average the various contributions are 
weighted with their quantum efficiency and that the average decay time < t > differs from 
p-r- (see Section E) . Indeed, for the data in Fig. 0](b) Eq. [21] yielded an average decay time 
of 31.1 ns , strongly different from the jf- — value of 10.4 ns, in contrast to the result (Eq. 
[5]) for single-exponential decay. 

It is important to note that Eq. [21] is the average decay time < t > corresponding to the 
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the average time given by Eq. 



decay curve given by Eq. [20] In, for instance, Refs. 
has erroneously been associated to fluorescence decay described by Eq. [TH] 

In contrast to the single-exponential model, the reliability function and the probability 
density function of a stretched-exponential do not have the same form (see Fig. [2]); the 
probability density function contains a time-dependent pre-factor. Therefore, the relation 
between the reliability function and the probability density function (Eq. [1]) has important 
consequences. For a /3-value of 0.5 the average decay time of the reliability function (Eq. 
T9l) and of the probability density function (Eq. [201 differ by more than a factor of ten. 
Thus it is important to take into consideration whether Eq. [19] or Eq. [20] is used to 
describe the experimental photoluminescence decay curve. This is important since in many 
reports^&i^L^ 1 ^, the luminescence decay curve in modelled with the time dependence 
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of Eq. [T9] We remark that while Eq. [19] can be used to account for the deviation from 
single-exponential decay, it does not represent the true Kohlrausch function, but is simply an 
alternative model. We argue that using the Kohlrausch stretched-exponential as a reliability 
function to model the fraction ^o) implies that the proper probability density function, i.e., 
Eq. [20], must be used to model a luminescence decay curve. Fig. IU(b) shows the modelling 
of experimental data with Eq. [20] with T str and (3 as adjustable parameters. The /3-value 
was 0.42 and T str was 96.2 /is -1 (p^- = 10.4 ns). Modelling with stretched-exponential is 
obviously more satisfactory than single-exponential, but here fails at long times, reflected 
by the high Xr _vame °f 60.7. 



III. CONCLUSIONS 



We have presented a statistical analysis of time-resolved spontaneous emission decay 
curves from ensembles of emitters, in particular colloidal quantum dots, with the aim to 
interpret ubiquitous non-single-exponential decay. Contrary to what is widely assumed, 
the density of excited emitters c(t) and the intensity in an emission decay curve (f(t) or 
g(t)) are not proportional, but the density is a time-integral of the intensity. The integral 
relation is crucial to correctly interpret non-single-exponential decay. We have derived the 
proper normalization for both a discrete, and a continuous distribution of rates, where 
every decay component is multiplied with its radiative decay rate. A central result of our 
paper is the derivation of the emission decay curve f(t) in case that both radiative and 
non-radiative decays are independently distributed (Eq. [T3|) . In this case, the well-known 
emission quantum efficiency can not be expressed by a single number anymore, but it is 
also distributed. We derive a practical description of non-single-exponential emission decay 
curves in terms of a distribution of total decay rates weighted with the radiative rates. 
Analyzing decay curves in terms of decay rate distributions opposes to the usual and widely 
reported analysis in terms of distributed lifetimes. We apply our analysis to recent examples 
of colloidal quantum dot emission in suspensions and in photonic crystals, and we find that 
this important class of emitters is well described by a log-normal distribution of decay 
rates with a narrow and a broad distribution, respectively. Finally, we briefly discuss the 
Kohlrausch stretched-exponential model; we deduce the average decay time and we find that 
its normalization is ill-defined for emitters with a realistic quantum efficiency of less than 
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100 % 
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